The momentum distribution of J/ip in B decays 
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The discrepancy between theory and data in the momentum distribution of slow J/t/j in B decays 
has been several times addressed as a puzzle. Using the most recent results on exclusive B decays 
into J/V' and heavy kaons or exotic mesons and reconsidering the non-relativistic-QCD calculation 
of the color octet fragmentation component, we show that an improvement in the comparison 
between data and theory can be obtained. There is still room for a better fit to data and this may 
imply that new exotic mesons of the XYZ kind have yet to be discovered. 

PACS; 14.40.Rt, 13.85.Ni, 13.25.Hw. 



In the first study of the inclusive decays of B mesons to J/V', the CLEO collaboration [T] observed that 
the J/ip momentum distribution dT/dp^ is well described at large by the sum of the exclusive modes 
B KJ/ip and B K*J/il} [2]. At lower momenta there are significant contributions from feed-down 
processes B — ^ XcX, Xc J/ip 7 and B — >■ ip'^: 4'' ~^ J l"^ 7r"'"7r~. After combining these contributions with the 
two-body modes, there naturally remains a shortfall across a large range of p^, having not counted contributions 
from higher kaon resonances and non-resonant multi-particle final states. 

Beneke et at [Ij proposed that the remaining part of spectrum could be explained by a large cc color octet 
contribution which feeds non-resonant multi-body final states; the calculation is performed within the framework 
of non-relativistic-QCD (NRQCD). They confront their results with the inclusive J/ip spectrum after having 
subtracted from it the B — >■ KJ/ip and B — K*J/%1] components. The agreement found with data is rather 
good using reasonable values for the physical parameters involved in the computation, essentially A^^^^^ and the 
Fermi momentum Pp of the b quark inside the B meson. 

Belle subsequently observed the decay B — >■ Ki{\21Q)J /il) with a branching fraction larger than B KJ /ij) 
and B — >■ K* J/tp. In an improved analysis of the inclusive J/ip spectrum, CLEO found that by summing all 
three of the two-body modes, a good fit to the inclusive spectra is obtained for p^, > 1.5 GeV 6J, which implies 
that the color octet contribution should be refitted. 




FIG. 1: The black-solid line represents the sum 
of all the contributions reported in Fig. |2j namely 
B — >■ ICJ/tp, B — ^ ICX and the octet contribution 
for Aqcj5 = 500 MeV and = 500 MeV, compared 
to the old theoretical prediction ^ (red-dashed line) 
computed as the sum of B ^ KJ/xp, B — >■ K*J/xp 
and the color octet component with Aq^j^ = 300 MeV 
and = 300 MeV. Going from the red-dashed line to 
the black-solid one the x^/DOF improves from 60/19 
to 22/19, having included also the theoretical errors. 
If one choses Aq^d = 800 MeV and p^ = 300 MeV 
the best fit further improves (x^/DOF = 19/19). 
Data points, black disks, are taken from BaBar [5]. 



FIG. 2: Different contributions to the inclusive spec- 
trum of B — >■ Jjip -\- All decays. The red-solid 
line accounts for the two-body decays of the type 
B — >■ K,,] /ip; the blue dot-dashed line shows the con- 
tributions from B — >■ KJ/tp and B — >■ K* J j-pi decays 
considered in [3] ; the green-dashed line represents the 
decays mediated by the exotic mesons B — >■ YZX — > 
ICJ /ip -\- light hadrons; the purple-dotted lines come 
from the non-resonant multi-particle final states, i.e., 
the color octet component, for three values of A^^^ = 
300, 500, 800 MeV (and p^ = 300 MeV) from right to 
left respectively, according to [J. Data points, black 
disks, are taken from BaBar [3]. 
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In a higher statistics analysis by BaBar [3], in which many feed-down modes from higher charmonia were 
directly measured and subtracted, the inclusive spectrum is confronted with the results of Beneke et al. , and it 
is noted that even after including color octet contributions, there remains an excess of events at low momenta: 
see Fig. [T] red-dashed curve in the region < < 0.8 GeV. This excess has provoked, since its observation, 
a variety of exotic interpretations. In [7] the discrepancy has been explained assuming that the J/^p recoils 
against a A — p strange baryonium state, while in [8l [9] the existence of intrinsic charm inside the _B-meson 
has been considered. Nevertheless the most explored possibility is the production of a strange hybrid meson 
Kh = sdg together with the J/ip in the B-meson decay. First proposed in [3], its contribution has been later 
quantified in [TDVfT^ . The BaBar analysis, however, does not include the large Ki{1270)J/tp mode, thus, for 
the reasons outlined above, it does not yield the full picture. 

More recently. Belle |13j analyzed the large i?+ — ^ J/V^ iiT+Tr+Tr" decay and performed an amplitude anal- 
ysis to determine the resonant structure of the K^tt^tt^ system, identifying several heavy kaon resonances 
/C = /i:i(1270),i^i(1400),A'*(1410),i4'2*(1430), 1^2(1600), X2(1770),i^2(1980). The two-body modes B ^ JCJ/^Ij 
further constrain the contribution from non-resonant (color octet) final states, and moreover, given that these 
newly identified kaons have large masses, they are expected to make significant contributions to the lower side 
of the momentum spectrum of J/ip- The red-solid line in Fig. [2] is our picture of all the above listed exclusive 
two-body contributions. Not including the color octet component, the agreement with data we found extends 
down to p.0 ~ 1.2 GeV. 

We also notice that XYZ exotic resonances are found to fill the low p^p bins as shown in Fig. [2j the green- 
dashed curve. Even if their weight in terms of branching ratio turns out to be modest with respect to what 
would be needed, their distribution peaks in the right region as commented in [14!. 

We thus reconsider the role of the NRQCD based model by Beneke et al. [4] with the aim of fully reconstructing 
the entire dT/dp^ distribution. In the following we discuss our computation of the black-solid curve in Fig. [l] 
which includes a new fit of the color octet component (see octet curves in Fig. [2]) and draw our conclusions. 
Essentially we notice that if new states of the XYZ kind will be found, since the fit turns out to be sensitive 
to this component, they would aid the color octet contribution to consistently fill the gap in the data. The 
observed relative small weight of XYZ in terms of branching ratios, as shown in Fig. [2] requires instead to force 
the NRQCD calculation in the high allowed region of values for Aq;,^ and p^. 

The heavy Kaons and XYZ contributions. Belle [13] measures the branching ratio 

Stot = B{B+ J/tA K+tt+'r-) = (71.6 ± 1 ± 6) x 10"^ (1) 

Looking at the invariant mass spectrum of X+tt+tt^, they are able to isolate the resonant contributions in 
which K^t:^'k~ originate from the decay of a heavy kaon K,j through some intermediate resonant state, namely 
7^, = Kp/uj,K*T:,K^{lA-iQ)'K,Kl{im)iT,Kfo^2, as in 

B{B+ ^ /C, J/^ ^ 7^, J/^ ^ J/V if+TT+TT-) = Btotfi (2) 

with fl the measured fractions. Since the interference between different heavy kaons is neglected, the sum of 
the measured fractions is larger than one. The neglected interference contributions are expected to be irrelevant 
when T]Cj/m]Cj << 1, so we rescale each fraction // proportionally to the width of the heavy kaon ICj 

where C is chosen in such a way that the sum over i and j of /| adds up to one. Now we observe that 

BiB+ JCjJ/ip n^J/i; -> J/^ K+tt+t:-) = Z, x B{B+ IC^J/^j) x B{ICj Ui) x B{n, Km:), (4) 

where li are isospin factors B{TZi Km) are all unity except for B{Ku! Ktttt) — 0.0153 and B{Kf2 — )• 
Knir) = 0.848. For some of the heavy kaons, namely i4:2(1600), iv:2(1770) and i^2(1980), the values of B{K.j ^• 
TZi) are not known experimentally and we thus extract a maximum value for them and in turn a minimum value 
for the relative B{B — ^ KjJ/tj)). The results obtained are summarized in Table IT] 



X{Kp) = l/2., X(i^*7r) =X(K'* 2(1430)77) = 4/9, X{Ku) = l, X{K fo) = X{K fi) = 2/3. 
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rriK, (GeV) 


Tk, (GeV) 


B{B+ JCjJ/^l)) X 10^ 


7^1(1270) 


1.270 


0.090 


144.0 ±29.3 


A'i(1400) 


1.403 


0.174 


25.1 ± 5.7 


A:*(1410) 


1.414 


0.232 


> 5.1 ±2.4 and < 11.8 ± 5.7 


A-2*(1430) 


1.430 


0.100 


40.2 ±24.0 


7^2(1600) 


1.605 


0.115 


> 8.4 ±2.9 


A'2(1770) 


1.773 


0.186 


> 4.4 ± 1.5 


^-2(1980) 


1.973 


0.373 


> 15.2 ±2.5 



TABLE I: Branching ratios for the decays B — ^ ICjJ/tp extracted from Belle [13] . 



The branching ratios in Table [T] are used to perform a Monte Carlo simulation of the decay chain, weight- 
ing vertices with the appropriate powers of momenta required by the decay partial wave ^. The sum of all 
contributions is shown by the red-solid curve of Fig. [2j 

We do the same for the exotic XYZ states. Even if almost all of these mesons decay into final states containing 
a charmonium state together with light hadrons, they cannot be easily interpreted as standard charmonia and 
thus are referred to as exotics |15| . Some of them have been observed in B decays produced together with the 
pseudoscalar kaon K. The relative branching ratios are reported in Table|Tl] In addition to these decays one can 





mx^ (GeV) 


Tx, (GeV) 


Final State 


B(B KXj KJ/'4> + hght hadrons) x 10^ 


X(3872) 


3.872 


0.003 


J/ip p — > J/ip 7r+7r~ 


0.72 ± 0.22 .16 


J/lp UJ 


0.6 ±0.3 [m 


F(3940) 


3.940 


0.087 


J/ij) 


3.70 ± 1.14 [H] 


F(4140) 


4.140 


0.012 




0.9 ±0.4 [19] 


F(4260) 


4.260 


0.095 




2.00 ±0.73 .18 



TABLE II: Measured branching ratios for the decay B — >■ KXj. Where data is available for both neutral and charged B 
we take the average. 



suppose that XYZ are produced with heavier kaons, when allowed by kinematics. Assuming that all the exotic 
particles X are spin 1 states (most of them do not yet have definitive quantum numbers) the transition matrix 
elements are given by {X{€,p)IC{q)\B{P)) — g e ■ q ioi a, spin kaon and by {X{e,p)IC{r],q)\B{P)) = g' e ■ rj for 
a spin 1 kaon"^. From a dimensional analysis [g] = M" and [g'] = M. Since we only have data on the branching 
ratios B{B — >■ KX) we need some hypothesis on the relationship between g and g' . We suppose that g' — A g, 
with A some mass scale. Using B{B KX{3872)) reported in Table |lT] and the upper limit |20j 

B{B -> K*X{3&72)) X B{X{3872) J/i; Tr+Tr") < 0.34 x 10"^ (5) 

we obtain A > 550 MeV. We thus decide to extract g from the known branching ratios reported in Table 11^ and 
fix g' = m^^j^-^^g^ where m]c{J — 1) is the mass of a generic spin 1 kaon. The green-dashed line in Fig. p shows 
the sum of all XYZ contributions along with all kaons allowed by kinematics. 

The color octet component. We reanalyzed the Beneke et al. 4J calculation to refit the color octet curve 
contribution to the J/ip spectrum. Their computation, based on the NRQCD factorization approach, assumes 
that the production of the cc pair in a color octet state, which has approximately the same probability as in 
color singlet, has a large kinematic effect on the momentum spectrum, due to the energy distribution of the soft 
gluons emitted by the {cc)s fragmentation into J/ip- 

The {cc)s pair can be produced, at order v"* in the non relativistic expansion, in three different L — S con- 
figurations: l^S'o, l^Pj and l^S'i. One can associate with each of these configurations a non-perturbative shape 
function for the energy and an invariant mass distribution of the radiated system (with a characteristic energy 
scale of rricv'^ w ^qcd)- The normalizations of the shape functions are related to the conventional NRQCD 



^ Since the heaviest kaons are spin 2 states the decay occurs in P-wave and thus events are naturally pushed towards higher values 

of in the allowed kinematic region. 
^ We cannot consider here the decay B — > K,X with a spin 2 kaon, since we do not have any experimental data to compare with. 
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matrix elements {On ), which describe the hadronization into a J/^p from a cc pair in an angular momentum 
and color state n. The {On^) can be determined by fits to J /ip production in a variety of processes. Nev- 
ertheless, enforcing these normalization conditions underestimates data, because the phenomenological values 
of the matrix elements are computed from integrated quantities in leading order calculations, while the shape 
functions contain higher order corrections in the velocity expansion. To overcome this difficulty one chooses 
to fix the absolute normalization by adjusting the sum of all contributions to data. The Fermi motion of the 
h quark inside the B meson is taken into account using the ACCMM model [H]. The tunable parameters are 
essentially A^^^^^ and the Fermi momentum . 

In this approach the color octet channel leads to non-resonant multi-body final states since the probability 
that the emitted soft gluons reassemble with the spectators (light quarks and gluons inside the B) to form a 
single hadron is assumed to be very small (factorization hypothesis). The two-body modes are due to the color 
singlet channel and must be considered separately. 

In [1] the comparison is made with data where the two-body modes with K and K* are subtracted. The 
agreement with the experimental spectrum, having fixed A^^j, = 300 MeV and = 300 MeV, is shown in 
Fig. [T] red-dashed curve. 

An observation is in order: possible interference effects between the octet final states, which at least contain 
two pions in addition to /C J/?/;, and the multi-body final states originated from the exclusive modes with heavier 
kaons decaying to /C4- pions and XYZ resonances can be safely neglected as Fyc^ /mx,j < 0.2 and Txj/mXj << 1. 
This would be a stronger assumption for hypothetical new broad resonances. 

Results and outlook. The black-solid curve in Fig. [l] is obtained as a sum of the standard two-body 
contributions from kaons (red-solid curve in Fig. [2]) with the XYZ contributions (green-dashed) plus one of the 
color octet components (purple-dotted). The latter are three curves obtained with A^^^^ = 300,500,800 MeV 
(and Pp — 300 MeV) from right to left respectively. To best fit data in Fig. [l] we choose the octet component 
with Aq^j, = 500 MeV and also need to push p^ up to p^, — 500 MeV. The values chosen for A^^j^ and Pp 
are critically on the high sides of the allowed ranges which are supposed to be A^^^^ g [200, 450] MeV [22] and 
Pf. £ [300, 450] MeV. Yet the black-solid curve in Fig. [l] represents a considerable improvement with respect 
to the old one (red-dashed). Relying on the validity of the NRQCD approach, our results seem to indicate 
that the addition of new resonances of the XYZ kind feeding the low p^ region would effectively improve the 
agreement with data. Indeed if the total branching ratio due to the XYZ turns out to be three times the 
currently observed value, a good description of data, namely x^/DOF = 28/19, would be obtained including a 
color octet component with A^^^^ = 300 MeV and p^ = 500 MeV. In this respect our results could be suggestive 
of the existence of a number of not yet discovered exotic mesons: new outcomes may arrive from Belle and 
LHCb. 
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